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ABSTRACT 

Atmospheric  density  fluctuations  induce  phase 
perturbations  that  degrade  the  spatial  coherence  of  a  laser 
beam  propagating  through  the  atmosphere.  These 
degradations  spread  the  laser  beam  and  alter  the  centroid 
and  intensity  profile  stochastically.  Turbulent  conditions 
arising  from  various  environmental  situations  and 
meteorological  phenomena  are  found  at  virtually  all  levels 
of  the  atmosphere.  A  substantial  fraction  of  the  optical 
turbulence  along  a  vertical  path  arises  from  the  surface 
heat  flux  within  the  first  100-200  meters  above  the  ground. 

This  thesis  seeks  to  measure  and  analyze  these 
turbulent  layers.  A  high  frequency  acoustic  echosounder 
was  developed  to  analyze  atmospheric  turbulence  within  the 
first  100-200  meters  above  the  ground  with  extremely  high 
accuracy.  The  echosounder  design  incorporated  a  25 
element,  square,  planar  array  housed  within  a  lead  and  foam 
insulated,  acoustic  shroud.  This  acoustic  echosounder  was 
employed  to  collect  real-time,  low  level  atmospheric  data 
at  two  sites.  This  data  should  prove  useful  in  determining 
whether  the  performance  of  laser  and  electro-optical 
systems  could  be  increased  by  raising  sensors  and 
transmitters  well  above  the  local  turbulent  layer  through 
the  use  of  towers  or  other  such  structures. 
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I.   INTRODUCTION 

A  coherent  laser  beam  propagating  through  the 
atmosphere  is  very  susceptible  to  numerous  turbulent- 
dependent  processes.  As  electromagnetic  waves  transit  the 
turbulent  regions,  atmospheric  irregularities  randomize  the 
amplitude  and  phase  of  the  wave.  A  high  resolution 
acoustic  profiler  was  designed  and  developed  to  quantify 
the  altitude  dependence  of  these  atmospheric  disturbances. 
Acoustic  energy  is  more  susceptible  to  these  irregularities 
than  electromagnetic  waves.  Historically  acoustic  sounders 
have  proven  to  be  the  best  means  of  providing  high 
resolution  atmospheric  turbulence  profiles  within  the  first 
500  to  1000  meters  above  the  ground.  Acoustic  echosounders 
are  frequently  used  to  detect  and  measure  atmospheric 
density  and  velocity  fluctuations  resulting  from  wind 
shears,  convection,  and  temperature  inversions. 

Presently  two  optical  parameters,  the  spatial  coherence 
length  ( ro )  and  the  isoplanatic  angle  (9o)  are  measures  of 
the  perturbation  of  electromagnetic  waves  propagating 
through  the  atmosphere,  and  are  measured  very  reliably  by 
optical  systems  developed  by  Walters  [Refs.  1  and  2]. 
Although  these  systems  measure  a  line  integral  of  the 
atmospheric  optical  turbulence  with  high  accuracy,  a  major 
drawback  of   each  system   is  that  no  provision   is  provided 


for  determing  the  height  of  the  atmospheric  disruptions. 
If  these  disturbances  are  found  to  exist  very  near  the 
surface,  it  may  be  possible  to  negate  their  effect  by 
elevating  sensors  and  transmitters  or  altering  the  surface 
heat  flux  to  suppress  the  generation  of  turbulence. 

This  thesis  deals  with  the  design,  construction  and 
implementation  of  a  high  frequency,  acoustic  echosounder 
which  will  accurately  analyze  the  atmospheric  density 
fluctuations  within  approximately  the  first  200  meters  of 
the  atmosphere.  As  this  project  is  a  product  of  the 
research  and  efforts  of  two  students,  the  work  was 
appropriately  divided.  My  particular  task  was  concentrated 
in  the  actual  design  and  hardware  development  of  the 
echosounder.  My  colleage,  LT .  Michael  Wroblewski  [Ref .  3], 
dedicated  his  efforts  toward  software  development  and  data 
display,  analysis  and  storage. 

The  echosounder  developed  and  discussed  herein  utilizes 
a  25  speaker,  square,  planar  antenna  array  enclosed  within 
an  acoustically  absorbent  shroud  comprised  of  lead  and  foam 
insulation  inside  a  sturdy  plywood  framework.  This  antenna 
array  was  employed  as  both  the  transmitter  and  receiver  of 
the  acoustic  signal.  A  computer  controlled  function 
generator  transmits  an  acoustic  pulse  through  the  antenna 
array.  Additional  electronic  components  incorporated 
within  the  overall  echosounder  system  are  used  to  amplify, 
digitize,  display,   analyze  and  ultimately  store  the  return 
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signal.  Over  time  the  received  signal  data  is  used  to 
produce  a  quantitative,  three  dimensional  map  of  the 
atmospheric  turbulence. 

Acoustic  echosounders  have  been  developed  and  in  use 
for  many  years  and  have  proven  to  be  a  valuable  probe  for 
analyzing  the  structure  and  dynamics  of  the  lower 
atmosphere  [Refs.  4  through  7].  Devices  similar  to  ours 
have  been  used  to  obtain  profiles  of  the  atmospheric 
density  and  temperature  fluctuations  [Refs.  8  and  9].  Our 
device  excels  by  utilizing  a  high  speed  HP217  computer  to 
control  and  monitor  the  echosounder.  This  enables  us  to 
obtain  real-time  data  collection  and  analysis  with  the 
ability  to  store  and  reproduce  the  atmospheric  profile 
plots  at  will.  This  information  should  help  quantify  the 
effects  of  low  level  turbulence  on  electro-optical  systems 
performance.  It  may  then  be  possible  to  negate  these 
turbulent  effects  by  simply  mounting  sensors  and 
transmitters  well  above  the  turbulent  layers  or  by  adaptive 
optics . 

This  thesis  discusses  the  development  of  a  high 
frequency,  5KHz ,  acoustic  echosounder,  the  criteria 
involved  in  the  actual  designs  of  both  the  acoustic  array 
and  the  acoustically  absorbent  enclosure,  and  the  selection 
and  function  of  all  the  associated  component  parts  of  the 
echosounder  system.  Finally  the  operational  success  of  the 
system  will  be  discussed  and  the  output  data  analyzed. 


II.   BACKGROUND 

Acoustic  echosounders  probe  the  atmosphere  by 
transmitting  a  pulse  of  acoustic  power  which  is 
subsequently  scattered  back  from  the  atmosphere  by 
temperature  and  velocity  inhomogeneities .  The  echosounder 
(echosonde)  equation,  often  referred  to  as  the  radar 
equation  in  meteorology,  is  used  to  determine  the 
backscattered  acoustic  power.  This  equation  is  summarized 
by  Neff  in  Reference  10  and  is  based  upon  the  work  of 
Tatarski  [Ref .  11]  and  Little  [Ref .  5]. 

Pr    =   Er[    PtEt     ][    exp(-2aR)    ][    ao(R,f)    ][    ^ct    ][    AGR- 2     ] 

where 

Pr  is  the  electrical  power  returned  from  a  range  R, 
Pt  is  the  electrical  power  transmitted  at  frequency  f , 
Er  is   the  efficiency  of  conversion  from  acoustic  power 

to  electrical  power  by  the  transducer, 

Et  is   the   efficiency   of   conversion   from  electrical 

power  to  acoustic  power  by  the  transducer, 

exp(-2aR)   is    the   round    trip   power   loss   due   to 

attenuation  where  a  is   the  average   attenuation  (meters-1 ) 

to  the  scattering  volume  at  the  range  R  (meters), 

ao(R,f)  is  the  scattering  cross  section  per  unit  volume 

at  a  distance  R  and  frequency  f, 
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c  is  the  local  speed  of  sound  (meters/second), 

t  is  the  pulse  length  (seconds), 

A  is  the  aperture  area  of  the  antenna  (meters2  ), 

R  is  the  range  (meters),  and 

G  is  the  effective-aperture  factor  of  the  antenna. 

Empirically  measuring  or  calculating  the  values  for  all 
other  terms,  one  can  use  this  equation  to  determine 
o~o(R,f),  the  scattering  cross  section  per  unit  volume;  that 
is,  the  fraction  of  incident  power  backscattered  per  unit 
distance  into  a  unit  solid  angle  at  a  frequency  f.  Based 
upon  experimental  results,  Tatarski  [Ref.  11]  expresses  the 
acoustic  backscatter  cross  section  per  unit  volume, 
ao(R.f),  by. 


Ct2 

ao(R,f)  -    0.0039  kW3   

To2 


where 

k   =   27r/A    is   the   incident   acoustic   wavenumber  at 

wavelength  A , 

To  is  the  local  mean  temperature  in  degrees  Kelvin,  and 
Ct2  is  the  temperature  structure  parameter. 

Combining  this  equation  with   the   echosonde   equation,  one 

obtains  a  volume-averaged  measure  of  Ct2. 


11  2  1  PR 

CT* To2    k-i/3 R2exp(2ocR) 

0.0039      ErEt  ct         AG        Pt 
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Hall  and  Wescott  [Ref.  12]  calculated  a  beam-shape 
compensation  factor  of  0.40  for  a  piston  source  antenna 
with  a  uniformly  illuminated  square  aperture.  This  value 
is  the  same  as  the  effective-aperture  factor,  G  and  can  be 
substituted  into  the  above  equation.  Approximating  the 
aperture  area  of  the  antenna  to  be  equivalent  to  25  times 
the  aperture  area  of  a  single  speaker  having  a  diameter  of 
7.620  centimeters,  we  get  a  value  of  0.1140  square  meters. 
Combining  these  values  with  the  numerical  constants  in  the 
above  equation,  we  can  simplify  the  equation  for  Ct2. 


1  1  PR 

Ct2    =    11245 To2    k"i/3 R2exp(2aR) 

Er  Et  ct        Pt 


The  efficiency  of  conversion  factors  for  our 
echosounder  were  determined  with  the  use  of  an  IVIE 
Electronics  IE-30A  Audio  Analyzer  Sound  Level  Meter,  a 
HP3314A  Function  Generator,  a  HP3561A  Dynamic  Signal 
Analyzer,  a  reference  speaker,  the  acoustic  array  and  the 
Naval  Postgraduate  School  anechoic  chamber.  The  function 
generator  was  used  to  transmit  a  1.0  volt,  5.0  KHz, 
sinusoidal  signal  through  the  reference  speaker.  The 
acoustic  array  was  employed  as  a  receiver  and  was  aligned 
along  the  centroid  of  the  main  lobe  of  the  reference 
speaker  signal.  With  the  acoustic  array  in  position,  the 
sound   level   meter   was   used  to  determine  the  total  power 
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available  for  reception  by  the  acoustic  array;  that  is,  the 
total  power  actually  ensonif ied  within  the  solid  angle  of 
the  acoustic  array.  Ten  independent  measurements  were  made 
with  the  sound  level  meter  and  then  averaged  to  yield  a 
power  value  of  73.7  dB.  Using  a  conversion  factor  and  the 
array  aperture  area  determined  above,  the  total  acoustical 
power,  Pa ,  available  at  the  acoustic  array  is  given  by  the 
equation  below. 

Pa  =  107.3  7  -  (  i  .  10-12  Watts  /  m2  )  -  (  0.11401  m2  ) 
or, 

Pa  =  2.67  •  10-6  Watts 

The  dynamic  signal  analyzer  was  then  used  to  measure  the 
actual  power  received  by  the  acoustic  array  and  converted 
into  electrical  power.  Four  independent  measurements  were 
recorded  and  averaged  to  yield  a  value  of  4.035  mVrms  .  The 
acoustic  array  impedance  was  previously  measured  and 
determined  to  be  12.1  ±0.4  Q.  These  values  were  used  in 
the  following  equation  to  determine  Pe ,  the  total 
electrical  power  converted  by  the  acoustic  array. 


(  4.035  •  10-3  Vrms   )2 
-  L34  .  io-6  Watts 

12.1  Q 


The   efficiency   of   conversion   factor   for    acoustic   to 
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electrical  power  is  the  ratio  of  Pa  to  Pe ,  or, 


1.34  •  10-6  Watts 

Er  = =  0.503 

2.67  •  10-6  Watts 


The  value  for  the  efficiency  of  conversion  for  electrical 
to  acoustic  power,  Et ,  is  not  readily  calculatable  with  the 
equipment  at  hand.  However,  based  upon  the  design  of  the 
speakers  used,  the  values  for  Er  and  Et  should  be 
approximately  equal.  Furthermore,  allowing  for 
experimental  error  especially  in  the  reading  of  the  sound 
level  meter,  the  best  accuracy  for  the  values  of  Er  and  Et 
is  0.5  ±0.1.  Substituting  conversion  factor  values  of  0.5 
into  the  above  equation,  we  can  further  simplify  the 
equation  for  Ct2 . 


To2        PR 
Ct2     =    44980    k-i/3 R2    exp(2aR) 

CT     PT 


The  above  equation  was  incorporated  into  a  computer 
program  that  provides  the  system  control,  data  acquisition 
and  data  processing  techniques  [Ref .  3].  The  reduced  data 
was  then  used  to  measure  the  temperature  structure 
parameter  as  a  function  of  time  of  day  and  altitude  for 
various  sites. 
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III.   SYSTEM  DESIGN  AND  EQUIPMENT  DEVELOPMENT 

A .   SPEAKERS 

In  the  design  of  our  echosounder,  it  was  determined 
that  a  rapid  decay  time  was  required  to  obtain  accurate 
short  range  information.  Weight  restrictions  involved  with 
equipment  transportation  plus  high  efficiency  in  the  2KHz 
to  5KH2  range  led  to  our  decision  to  use  piezo  ceramic 
speakers.  The  Motorola  KSN  1005A  speaker  was  selected 
based  on  these  requirements  and  the  specifications  charted 
in  the  Motorola  catalog  [Ref.  13]  and  reproduced  below  in 
Figures  1  through  3. 

< 
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Fig.   1.   Speaker  Ratings 
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Typical  Frequency  Response 
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Fig.   2.   Speaker  Frequency  Response 

Dimensions:  KSN  1005 A,  KSN  1003A 
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Speaker  Dimensions 
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These  charts  indicate  that  a  maximum  response  for  our 
speakers  occurs  at  a  resonant  frequency  of  5000  Hertz. 
This  frequency  was  used  as  the  baseline  from  which  all  our 
measurements  were  made. 

Using  the  speakers  in  an  anechoic  chamber,  the  average 
e- !  voltage  decay  time  was  measured  to  be  approximately  900 
usee  (Fig.  4).  This  decay  time  translates  into  a  sound 
propagation  distance  of  just  over  15.0  centimeters  (at  STP) 
from  the  speakers.  Considering  our  requirements,  this 
speaker  was  ideally  suited  to  serve  our  purpose. 
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Fig.   4.   Speaker  Decay  Time  Trace 
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B.   ACOUSTIC  ARRAY 

Our  next  consideration  was  the  echosounder  array 
pattern.  Ideally,  the  acoustic  sources  should  be  placed 
exactly  one  half  wavelength  apart.  At  a  frequency  of  5000 
Hertz,  this  would  require  spacings  of  3.4  centimeters  (at 
STP)  which  is  physically  impossible  for  the  speakers  we 
have  chosen.  The  closest  possible  spacing  was  7.62 
centimeters  between  sources  after  shaving  off  the  flange  of 
the  horn  (Fig.  3) . 

From  Kinsler,  et  al .  [Ref.  14:  p.  195],  the  equation 
for  the  directionality  factor  of  a  simple  line  array  is 
derived  as : 


H(N,8)  = 


N 

sin  (  kd  sin  9  ) 

1  2 

N  1 

sin  (  kd  sin  8  ) 

2 


where 

k  is  the  wave  number  ( 2n/  A.  ), 

d  is  the  distance  between  sources, 

N  is  the  number  of  sources,  and 

9  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest. 

However,  this  equation  assumes  simple  point  sources 
which  does  not  accurately  describe  the  speakers  we  have 
chosen.    It  was   necessary  to   couple  this  equation  to  the 
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directionality   factor   for   a   piston  source  which  is  also 
identified  in  Kinsler,  et  al .  [Ref .  14:  p.  108]. 


D(9)  = 


2  Ji  (  ka  sin  9  ) 
ka  sin  9 


where 

k  is  the  wave  number, 

a  is  the  radius  of  the  piston  source, 

9  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest,  and 

Ji  is  a  first  order  Bessel  function. 

These  two  equations  were  combined  to  produce  the 
equation  of  directionality  for  a  linear  array  of  piston 
sources,  L(N,9),  by  simple  multiplication. 


L(N,9)  =  H(N,9)  •  D(9) 


This  equation  was  incorporated  into  a  computer  program 
[Appendix  A]  and  used  to  generate  linear  array  patterns  for 
varying  numbers  of  piston  sources  (Figs.  5  through  9). 
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PATTERN  FOR   3   ELEMENTS 


PATTERN  TOR   4   ELEMENTS 


Fig.  5.  Three  Element  Array       Fig.  6.  Four  Element  Array 

PRTTERN  FOR   5   ELEMENTS 


Fig.   7.  Five  Element  Array 
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PATTERN  FOR  7      ELEMENTS 
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Fig.  8.  Six  Element  Array 


Fig.  9.  Seven  Element  Array 
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It  was  concluded  from  these  plots  of  the  linear  array  beam 
patterns  for  varying  numbers  of  piston  sources  that  five 
linear  elements  would  provide  the  best  combination  of 
forward  directionality,  sidelobe  suppression,  physical 
size,  and  a  relatively  low  cost.  Then,  in  order  to  enhance 
both  array  efficiency  and  symmetry,  we  settled  on  a  five  by 
five  element  array  design  with  7.62  centimeter  spacing 
between  speakers  in  both  the  vertical  and  horizontal 
directions . 

After  we  verified  the  manufacturer's  polarity 
designation  for  35  speakers  and  inspected  each  of  these 
speakers  for  physical  defects  and  loose  wiring  connections, 
we  obtained  Lissajous  plots  for  each  individual  speaker. 
Figure  10  represents  a  typical  Lissajous  plot  for  one  of 
the  speakers  incorporated  into  the  five  by  five  element 
array.  Based  upon  the  speaker's  output  to  input  voltage 
ratios  as  illustrated  by  these  Lissajous  plots,  we  were 
able  to  rank  all  of  our  speakers  by  signal  efficiency.  It 
was  based  upon  this  criteria  that  we  selected  the  25  most 
efficient  speakers  for  use  in  the  five  by  five  element 
array,  placing  the  best  speakers  at  the  center  and 
subsequent  ranking  speakers  further  toward  the  sides  and 
corners . 


21 


t 

Test 
Outpu 

speaker 

A 

-^ 

rence  Speaker  1 

> 

■ 

Fig.   10.   Lissajous  Plot 

The  25  selected  speakers  were  mounted  in  a  five  by  five 
planar  array  on  a  balsa  wood  insulated  bilayered  sheet 
metal  board.  After  wiring  all  the  speakers  in  series,  we 
surrounded  all  the  electrical  connections  and  speaker  backs 
with  two  3.0  centimeter  layers  of  foam  insulation 
sandwiching  a  1.0  millimeter  lead  sheet.  Then  the  entire 
array  mounting  was  enclosed  in  a  44  by  44  by  5  centimeter 
sheet  metal  box.  This  design  was  chosen  to  suppress 
virtually  all  acoustic  energy  propagating  out  the  rear 
hemisphere  of  the  array,  while  shielding  the  array  from  any 
external  electrical  interference  (Figs.  11  and  12). 
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Fig.   11.   Array  Photograph 
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Fig.   12.   Array  Photograph 

C.   ENCLOSURE 

Acoustic  echosounding  has   proven   to   be   an  extremely 
useful   technique   for   probing   and   analyzing   the   lower 
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atmosphere.  In  order  to  utilize  the  acoustic  waves 
transmitted  and  later  received  by  this  remote  sensing 
method  efficiently,  it  is  essential  to  have  an  efficient 
antenna  with  highly  directive  beams  and  strongly  suppressed 
sidelobes.  Antenna  design  becomes  increasingly  more 
important  in  a  noisy  environment  where  noise  pollution 
within  the  sidelobes  may  dominate  the  desired  signal  within 
the  main  lobe.  Hall  and  Wescott  [Ref .  12]  showed  that 
sidelobe  suppression  improved  with  higher  frequencies. 
Their  studies  showed  that  the  measured  90  degree  sidelobe 
suppression  ranged  from  38  dB  at  1  KHz  to  50  dB  at  5  KHz. 
Furthermore,  any  significant  improvement  in  sidelobe 
suppression  could  only  be  obtained  by  surrounding  the 
antenna  with  an  acoustic  energy  absorbing  cuff  or  shroud. 

In  an  effort  to  maximize  our  antenna  main  lobe  to 
sidelobe  power  ratio,  we  intend  to  operate  at  5KHz . 
Additionally,  we  have  designed  an  acoustic  energy  absorbing 
enclosure.  Many  designs  were  considered  based  upon 
previous  research  in  the  field  of  echosounding  [Refs.  15 
through  18].  In  addition,  we  obtained  the  actual  acoustic 
beam  patterns  for  our  array  using  a  computer  program 
written  by  LCDR  Butler  [Ref.  19]  which  we  modified  for  our 
purposes.  This  modified  version  of  LCDR  Butler's  computer 
program  is  in  Appendix  B.  By  rotating  the  array  in  an 
anechoic  chamber,  we  were  able  to  produce  highly  accurate 
polar  plots  of  the  array  beam  patterns  (Figs.  13  and  14). 
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5X5   ACOUSTIC   ARRAY   BEAM   PATTERN 


^Imike 


R=3.60  M         MIKE  VT=  8. 2465       INPUT-5. 0  V         FREQ=5.0KHZ 


Fig.       13.       Polar   Plot   of    5   X    5   Array   Beam   Pattern 
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DIAGONAL  '5   X   5  ACOUSTIC   ARRAY   BEAM   PATTERN 


R=3.  60  M         MIKE  VT=»  3.9497        INPUT-2.  5  V         FREQ=5.  0  KHZ 


Fig.       14.       Polar   Plot   of    5   X   5   Array   Beam   Pattern 

at    a   45°    Aspect 
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A  comparison  of  the  theoretical  beam  pattern  for  the 
five  element,  piston  source  array  in  Figure  7  with  the 
experimentally  measured  polar  plot  of  the  5  X  5  acoustic 
array  in  Figure  13  and  the  associated  coordinates 
corresponding  to  the  individual  data  points  illustrate 
significant  agreement.  As  predicted  by  our  computer  model 
[Appendix  A] ,  the  actual  acoustic  beam  pattern  is  comprised 
of  nine  distinct  lobes  in  the  forward  hemisphere.  Both  the 
positions  of  these  lobes  and  their  relative  amplitudes  are 
within  2.0  percent  of  their  predicted  values.  The  computer 
modeled  main  lobe  full  width  angle  of  19.5  degrees  very 
closely  predicted  the  experimentally  measured  main  lobe 
full  width  angle  of  19.0  degrees.  Reflected,  scattered, 
diffracted  and  electronic  noise  all  contribute  to  the 
clutter  near  the  origin  of  Figure  13.  This  is  the  primary 
reason  why  the  minima  don't  actually  go  to  zero  as 
predicted  by  our  model.  Additional  reasons  include 
mechanical  noise  in  the  device  used  to  rotate  the  acoustic 
array  and  vibrational  noise  in  the  array  itself.  The 
apparent  lobe  in  the  upper  left  hand  quadrant  is  an  anomaly 
seen  in  all  polar  plots  herein,  and  corresponds  to  acoustic 
energy  reflected  off  the  door  and  window  of  the  anechoic 
chamber.  Figure  14  represents  the  acoustic  array  mounted 
at  a  450  aspect.  A  computer  model  for  this  orientation  is 
a  very  complicated  task  and  was  not  performed  at  this  time. 
However,   this   polar   plot   of   the   acoustic   array   beam 
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pattern,  and  in  particular,  the  appreciable  sidelobe 
suppression  strongly  support  the  need  for  further  array 
design  studies. 

Since  it  was  our  aim  to  suppress  all  sidelobes  and 
utilize  the  main  lobe,  we  chose  not  to  taper  our  enclosure 
as  most  previous  researchers  had.  Rather  we  designed  the 
enclosure  based  upon  the  dimensions  of  the  array  itself  and 
the  acoustic  beams  it  generated.  The  height  of  the 
enclosure  is  160  cm.  Allowing  for  the  base  foam  and  lead 
insulation,  the  array  housing  and  the  base  framework 
itself,  the  actual  array  piston  sources  are  positioned  143 
cm  from  the  mouth  of  the  enclosure.  The  length  and  width 
of  the  enclosure  are  both  65  cm.  Again  allowing  for 
construction  materials,  the  inner  dimensions  are  both  48 
cm.  Hence,  the  position  of  the  center  speaker  is  24  cm 
from  any  of  the  four  sides.  These  dimensions  produce  the 
unobstructed  19.0  degree  path  into  which  the  main  lobe  is 
transmitted . 

2  •  tan-i  (  24  cm  /  143  cm  )  =  19.050  *  19. 0© 

Plywood  was  used  for  the  construction  of  the  enclosure 
and  provided  not  only  a  rigid,  inexpensive  framework,  but 
also  provided  an  impedance  mismatch  which  helped  to 
attenuate  external  noise  interference.  Anticipating  all 
kinds  of  weather  conditions  during  data  collection,  the 
plywood  enclosure  was  first  waterproofed  with  four  coats  of 
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marine   varnish.     Grooved   joints,   caulking   and  weather 
stripping  were  also  design  considerations. 

A  millimeter  layer  of  lead  can  suppress  an  acoustic 
signal  by  reflection  as  much  as  40  dB  (Figs.  15  and  16). 
About  7.0  centimeters  of  convoluted  foam  can  suppress  a 
signal  another  3  to  4  dB  (Figs.  15  and  17).  The  primary 
purpose  of  the  convoluted  foam  is  to  attenuate  the 
sidelobes  which  repeatedly  reflect  off  the  enclosure 
interior  at  glancing  incident  angles.  Together  the  lead  and 
convoluted  foam  make  an  extremely  efficient  absorbing 
material  for  use  in  our  enclosure. 
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Fig.   16.   Signal  Suppression  by  Lead 
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Fig.   17.   Signal  Suppression  by  Foam 
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Lead-lined  absorbing  foam  was  glued  to  all  inner 
surfaces  of  the  enclosure  with  two  1.0  millimeter  layers  of 
lead  overlapping  at  all  corners.  Strong  aluminum  brackets 
were  used  to  connect  the  four  side  panels  to  each  other  as 
well  as  to  the  enclosure  base  (Figs.  18  through  20). 
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Fig.   18.   Enclosure  Photograph,  Fully  Assembled 
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Fig.   19.   Enclosure  Photograph,  Interior 
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Fig.   20.   Enclosure  Photograph,  Base  with  Array 


34 


D.   ADDITIONAL  HARDWARE 

In  addition  to  the  acoustic  array  and  enclosure 
discussed  in  the  previous  sections,  there  are  primarily  six 
more  pieces  of  equipment  essential  to  the  function  of  the 
echosounder.  Below  is  a  brief  description  of  these 
additional  components . 

1 .  HP2Q0  Series  Computer 

The  Hewlett-Packard  (HP)  200  Series  Computer 
includes  a  20  megabyte  hard  disk  drive,  a  floppy  disk  drive 
and  an  associated  printer  and  monitor.  The  HP200  Series 
Computer  is  the  central  control  component  for  the  entire 
echosounder  arrangement.  The  computer  used  is  an  HP217 
programmed  in  Basic  3.0  and  equipped  with  an  Infotek  BC203 
Basic  Compiler  and  an  Infotek  FP210  Floating  Point 
Accelerator  to  enhance  the  speed  of  program  execution.  The 
computer  is  used  to  maintain  timing  and  process  control ,  to 
define  all  parameters  for  the  HP3314A  Function  Generator  as 
well  as  to  execute  the  trigger  command  that  produces  the 
echosounder  transmitted  acoustic  pulse.  The  computer  also 
processes  the  data  from  the  acoustic  array  via  the  pre- 
amplifier, bandpass  filter  and  analog  to  digital  converter. 
The  computer  then  conducts  data  reduction  routines  to 
produce  and  display  a  high  resolution  atmospheric  profile. 

2.  HP3314A  Function  Generator 

The  Hewlett-Packard  (HP)   3314A   Function  Generator 
is   a   multimode   function   generator   capable  of  providing 
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sine,  square  and  triangular  wave  functions  as  well  as  any- 
desired  waveform  ranging  in  frequency  from  0.001  Hertz  to 
19.999  Megahertz.  The  HP3314A  Function  Generator  was  used 
to  supply  a  pulse  of  an  integer  number  of  sinusoidal  cycles 
of  constant  amplitude  to  the  QSC  Model  1700  Audio 
Amplifier.  Typically  a  pulse  of  100  sinusoidal  cycles  was 
used.  A  constant  frequency  setting  of  5000  Hertz  was  used 
for  all  data  runs. 

3.   QSC  Model  1700  Audio  Amplifier 

The  QSC  Model  1700  Audio  Amplifier  is  a  high  power 
amplifier  which  can  supply  350  watts  over  an  8  ohm  load. 
This  amplifier  was  used  to  boost  the  output  voltage  of  the 
function  generator  by  a  factor  of  20  from  1 . 5  volts  to  30 
volts.  The  power  supplied  to  each  speaker  in  the  array 
during  operation  was  then  37.1  watts  as  indicated  in  the 
equation  below. 


(   Voltagerms  )2        (  30  )2 

Power  = = =  37.1  Watts 

Impedance  2  •  12.1 


4 .   Pre -amplifier 

The  pre-amplif ier  was  designed  and  constructed  by 
Walters  [Ref .  20]  to  supply  a  gain  of  1000  to  the  returned 
signal  and  isolate  the  data  aquisition  system  from  the 
transmitted  pulse.  The  pre-amplif ier  may  be  thought  of  as 
a  safety  and  switching  mechanism  for  the  system.  An  LT 
1007  and  an  LT  1037  Operational  Amplifier  were  selected  for 
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use  in  the  pre-amplif ier  based  upon  their  low  noise 
properties  which  were  evaluated  by  physically  incorporating 
them  into  the  device  and  measuring  their  noise  and  gain 
characteristics  with  a  HP3561A  Dynamic  Signal  Analyzer. 
The  IN  4000  Series  Rectifiers  isolate  the  power  amplifier 
from  the  receiver  between  pulses  and  limit  the  voltage 
applied  to  the  operational  amplifiers.  The  electrical 
diagram  of  this  device  is  outlined  in  Figure  21. 
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Fig.   21.   Schematic  of  Pre-amplif ier 
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5 .   Rockland  Wavetek  Model  852  Bandpass  Filter 

The  Rockland  Wavetek  Model  852  Filter  operates  as  a 
48  dB  per  octave  bandpass  filter  that  suppresses  broadband 
noise  in  the  system.  High  and  low  bandpass  settings  of 
5500  Hertz  and  4500  Hertz  respectively  were  used  with  a 
roll  off  of  48  dB  per  octave.  The  filter  response  at  these 
settings  is  illustrated  in  Figure  22. 
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Fig.   22.   Bandpass  Filter  Response 


The  3dB  (half  power)  bandwidth  was  approximately  1  KHz 
wide.  Although  it  would  have  been  desirable  to  reduce  the 
bandpass  interval  to  between  50  and  100  Hertz  without 
impairing  signal  response,  the   Rockland  Wavetek   Model  852 
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Bandpass   Filter    was   the   best   filter   available.     It 
performed  adequately  despite  the  large  1  KHz  bandpass. 
6 .   Infotek  AD200  Analog  to  Digital  Converter 

An  Infotek  AD200  12  bit  Analog  to  Digital  Converter 
was  used  to  sample  the  amplified  signal  voltage  and 
digitize  it  for  the  computer  system.  A  typical  sampling 
frequency  used  was  20000  Hertz. 

The  arrangement   of  the   above  equipment  is  illustrated 
in  Figure  23. 
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Fig.   23.   Echosounder  Equipment  Arrangement 


39 


E .   SOFTWARE 

The  software  and  the  HP217  Computer  were  responsible 
for  controlling  and  monitoring  every  operational  aspect  of 
the  system.  LT  Wroblewski  [Ref .  3]  provides  an  in-depth 
discussion  of  the  echosounder  software  which  is  also 
summarized  in  Appendix  C. 
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IV.   DATA  ANALYSIS 

A.   ECHOSOUNDER  PERFORMANCE 

An  analysis  of  the  acoustic  echosounder  output  was 
conducted  to  determine  the  validity  of  previously 
determined  echosounder  parameters.  The  e~ i  decay  time 
constant  for  the  applied  voltage  of  a  single  speaker  was 
previously  calculated  to  be  900  us  which  corresponds  to  a 
distance  of  15.3  cm.  The  e~  i  decay  time  constant  for  the 
applied  voltage  of  the  entire  acoustic  array  was  determined 
by  the  same  experimental  procedure  to  be  1060  us  which 
corresponds  to  a  distance  of  18.0  cm.  Although  the  25 
speakers  in  the  acoustic  array  transmitting  in  unison 
should  have  approximately  the  same  e- 1  decay  time  constant 
as  the  single  speaker,  it  is  not  surprising  that  the  array 
e~ i  decay  time  constant  is  larger.  Using  ten  time 
intervals,  an  expectation  of  99.99  percent  decay  of  the 
transmitted  signal  should  occur  within  a  distance  of  180.0 
cm  or  just  under  2.0  meters.  Reviewing  typical  echosounder 
traces  and  Ct 2  plots  [Appendix  D] ,  it  is  evident  that  such 
a  short  recovery  time  is  actually  quite  accurate.  The 
recovery  distance  is  typically  found  to  be  on  the  order  of 
9  meters  or  53  ms .  This  "blind  zone"  is  not  uncommon  in 
most  radar  systems  and  is  in  part  a  return  signal  from 
ground   clutter   surrounding   the   acoustic   array   and  its 
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enclosure  immediately  following  pulse  transmission.  An 
analog  to  digital  intensity  factor  of  4000  is  assigned  to 
the  heavily  shaded  region  of  the  recovery  zone.  Since  the 
computer  software  has  generated  this  value  as  the  product 
of  the  actual  analog  to  digital  "count"  value  and  the  range 
from  the  acoustic  echosounder,  this  value  of  4000  must  be 
divided  by  the  range  at  the  recovery  zone  boundary  to  give 
us  an  actual  analog  to  digital  "count"  value.  This  result 
is  then  multiplied  by  an  experimentally  measured  conversion 
factor  of  3.34  •  10-8  voltsr ms /count .  This  result  may  then 
be  used  in  the  determination  of  the  total  electrical  power 
received  and  converted  by  the  acoustic  array.  The  value 
for  the  total  electrical  power  received  by  the  acoustic 
array  is  calculated  in  the  equations  below. 


Voltagerms  =  (3.34  •  10-3v/count)  ( 4000m«counts )  (9m)-i 
Voltagerms  =  1.48  •  10~5  volts 

Power  =  (Voltagerms)2  •  (Array  Impedance)-1 
Power  =  (1.48  •  10-8  volts)2  •  (12.1  Q)-i 
Power  =  1.82  •  10~ H  Watts 


Using  this   received  electrical   power  value   and  the  known 
transmitted   acoustic   power   value   of   37.1  Watts,  we  can 
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solve  for  the  e~  1  decay  time  constant  of  the  applied 
voltage.  The  product  of  the  total  acoustic  power  received 
and  the  efficiency  of  conversion  factor  (  Er  =0.50  ), 
yields  the  total  electrical  power  received.  Hence,  the 
total  acoustic  power  received  is  3.64  •  10~  i  i  Watts. 
Although  the  distance  of  9  meters  corresponds  to  a  time 
interval  of  53  ms ,  this  value  is  actually  inflated  by  the 
system  pulse  length  of  20  ms .  Compensating  for  this,  the 
actual  recovery  zone  time  interval  is  33  ms .  The  e- 2  decay 
time  constant  of  the  power  is  equal  to  the  e_1  decay  time 
constant  of  the  applied  voltage.  These  values  are 
calculated  below. 


3.64  •  10-n  Watts 
-  e-2t/x   or, 

37. 1  Watts 


-  2  •  33  ms 

x    — =  2  .  56  ms   , 

In  (6.37  •  10-12 ) 


This  e-i  decay  time  constant  of  the  applied  voltage 
corresponds  to  a  distance  of  43.5  cm  and  is  about  twice  the 
previously  determined  value. 

Consideration  was  also  given   to  the   construction  of  a 
hexagonal  planar   array  rather   than  a   square  planar  array 
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based  upon  the  astounding  sidelobe  suppression  for  the  five 
by  five  element,  square  array  at  a  45°  orientation  as 
illustrated  in  Figure  14.  To  investigate  the  potential  of 
such  a  close  packed  design,  two  small  seven  element 
hexagonal  arrays  were  constructed.  The  first  array  design 
utilized  a  close  packed  spacing  of  7.62  cm,  while  the 
second  array  design  employed  a  3/2  wavelength  or  10.2  cm 
spacing.  Polar  plots  of  the  beam  patterns  for  both  these 
hexagonal  array  designs  were  made  by  rotating  the  arrays  in 
an  anechoic  chamber  with  the  main  diagonal  both  horizontal 
and  vertical.  In  both  cases  the  first  sidelobe  suppression 
for  both  close  packed  orientations  was  down  by  nearly  15 
dB,  while  the  first  sidelobe  suppression  for  both  3/2 
wavelength  spacing  orientations  was  only  down  by  about  7 
dB.  Furthermore,  the  ratios  of  the  size  of  the  main  lobe 
to  those  of  the  sidelobes  is  much  greater  for  the  close 
packed  design  than  they  are  for  the  3/2  wavelength  spacing 
arrangement  for  both  the  vertical  and  horizontal 
orientations.  The  polar  plots  of  the  array  beam  patterns 
for  both  close  packed  and  3/2  wavelength  spacing  hexagonal 
designs  are  illustrated  in  Figures  24  through  27.  These 
figures  strongly  warrant  further  study  of  the  hexagonal 
close  packed  array  design. 
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HORIZONTAL  DIAGONAL  CLOSE  PACKED.  7  SPEAKER  HEXAGONAL  ARRAY 


MIKE 


R-3.  60  M    MIKE  VT-  5. 6364   INPUT-5. 0  V    FREQ-5. 0  KHZ 


Fig.   24.   Polar  Plot  of  Close  Packed  Hexagonal  Array 

(Horizontal  Diagonal) 
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VERTICAL  DIAGONAL.  CLOSE  PACKED.  7  SPEAKER  HEXAGONAL  ARRAY 


MIKE 


R=3.60  M         MIKE  VT-  5.  1559       INPUT-5.0V         FREQ=5.0KHZ 


Fig.       25.       Polar  Plot   of    Close   Packed   Hexagonal   Array 

(Vertical    Diagonal) 
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HORIZONTAL  DIAGONAL    3/2  WAVELENGTH  SPACING.    7  SPEAKER  HEXAGONAL  ARRAY 


MIKE 


R-3.  60  M         MIKE  VT-  5. 18       INPUT-5. 0  V         FREQ-5. 0  KHZ 


Fig.       26.       Polar   Plot    of    3/2   Wavelength  Hexagonal   Array 

(Horizontal   Diagonal) 
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VERTICAL  DIAGONAL.  3/2  WAVELENGTH  SPACING.  7  SPEAKER  HEXAGONAL  ARRAY 


MIKE 


R=3.60  M    MIKE  VT-  6.0215   INPUT»5.0  V    FREQ=5.0  KHZ 


Fig.   27.   Polar  Plot  of  3/2  Wavelength  Hexagonal  Array 

(Vertical  Diagonal) 
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B.   SITE  EVALUATION 

Echosounder  data  was  collected  at  two  different 
locations.  The  primary  data  collection  site  was  the  upper 
roof  of  Spanagel  Hall  at  the  Naval  Postgraduate  School, 
Monterey,  California.  This  site  was  chosen  simply  for 
convenience  of  data  collection.  Data  gathered  at  this 
location  is  believed  to  represent  the  California  coast 
during  the  spring  near  sea  level.  The  second  data 
collection  site  chosen  was  in  the  vicinity  of  the  24  inch 
telescope  at  Lick  Observatory,  San  Jose,  California.  This 
site  is  located  atop  Mt .  Hamilton  at  an  altitude  of 
approximately  5700  feet  and  nearly  20  miles  inland  from  the 
coast . 

These  two  data  collection  sites  represent  areas  of 
differing  atmospheric  air  pressures,  water  vapor  pressures, 
local  temperature  ranges,  and  local  wind  velocity  ranges. 
These  characteristics  all  play  important  roles  in  effecting 
the  local  atmospheric  turbulent  conditions  and  thereby  the 
atmospheric  structure  parameter,  Ct 2 . 

In  addition  to  collecting  echosounder  data  at  Lick 
Observatory,  simultaneous  measurements  of  the  isoplanatic 
angle  ( 9o )  and  spatial  coherence  length  ( ro )  were  made  with 
systems  developed  by  Walters  [Refs.  1  and  2].  A  basic 
knowledge  of  these  two  systems  is  necessary  to  understand 
the  correlation   procedures   made.   The   isoplanatic   angle 
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(8o )  is  primarily  an  upper  atmospheric  measurement  which 
indicates  atmospheric  disruptions  at  a  range  of  two  to  15 
kilometers  above  the  Earth's  surface.  On  the  otherhand, 
the  spatial  coherence  length  (ro)  is  a  uniform  measure  of 
the  effects  of  the  entire  atmospheric  blanket  on  coherent 
light  transmission.  A  close  comparison  of  the  acoustic 
echosounder  data,  the  isoplanatic  angle  (80)  measurements 
and  the  spatial  coherence  length  ( ro )  measurements  should 
give  us  an  indication  of  the  turbulent  contributions  from 
the  lower  and  upper  troposphere  as  well  as  the  stratosphere 
above . 

Based  upon  isoplanatic  angle  ( 9o )  and  spatial  coherence 
length  ( ro )  measurements  at  Mt .  Wilson  near  Los  Angeles, 
California,  a  strong  correlation  between  the  isoplanatic 
angle  (9o)  and  spatial  coherence  length  ( ro )  measurements 
occurs  if  the  low  altitude  boundary  layer  contribution  is 
sufficiently  small.  This  strong  correlation  helps  to 
reinforce  the  overall  description  of  the  atmosphere  at  the 
time  of  data  collection.  Figures  28  and  29  graphically 
illustrate  the  atmospheric  measurements  made  at  Mt .  Wilson 
on  2  April  1987  by  Walters.  The  strong  correlation  between 
the  isoplanatic  angle  (9o)  and  the  spatial  coherence  length 
( ro )  is  especially  evident  between  the  hours  of  0700  and 
1300  universal  time.  The  close  tracking  of  these  two 
measurements   during  this   time  interval   indicate  that  the 
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Fig.   28.   Isoplanatic  Angle  Measurements,  Mt .  Wilson 
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Fig.  29.  Spatial  Coherence  Length  Measurements,  Mt .  Wilson 
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upper  atmospheric  conditions,  as  measured  by  the 
isoplanatic  angle  (60),  are  dominating  the  entire 
atmospheric  profile  as  measured  by  the  spatial  coherence 
length  ( ro ) .  Unfortunately,  the  existence  or  non-existence 
of  any  turbulent  surface  effects  cannot  be  ascertained  by 
the  employment  of  the  above  two  systems  alone.  However, 
employment  of  these  two  systems  together  with  the  acoustic 
echosounder  should  enable  us  to  produce  a  complete 
atmospheric  profile  with  strong  correlation  between  all 
three  atmospheric  measurements. 

On  9  and  10  April  1987,  all  three  systems  were  operated 
at  Lick  Observatory.  Again  good  correlation  between  the 
isoplanatic  angle  (60 )  and  the  spatial  coherence  length 
( ro )  measurements  was  noted.  In  addition,  strong 
correlation  between  the  spatial  coherence  length  ( ro ) 
measurements  and  the  acoustic  echosounder  data  was 
conclusively  present.  A  comparison  of  the  isoplanatic 
angle  (80)  and  the  spatial  coherence  length  ( ro ) 
measurements  in  Figures  30  and  31  shows  good  correlation 
between  the  two  atmospheric  parameters  especially  during 
the  0930  to  1130  time  interval  on  9  April  1987.  However, 
during  subsequent  hours  the  isoplanatic  angle  ( bo ) 
measurements  remain  high  (approximately  12  urad)  indicating 
relatively  calm  turbulent  conditions  in  the  upper 
atmosphere,  while  the  spatial   coherence  length  ( ro )  values 
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drop  sharply  after  1100  universal  time  indicating  dominant 
and  increasing  lower  atmospheric  turbulence.  This  increase 
in  the  low  level  turbulence  should  be  evident  in  the 
echosounder  data  commencing  around  1200  univers-al  time 
(0400  local  standard  time)  on  9  April  1987.  A  comparison 
of  the  echosounder  data  in  Figures  32  through  34 
illustrates  this  increase  in  the  local  surface  turbulence. 


3  Rpr  i  1  1937  Lick 


3      10     11     13 
Time  (hrs  UT) 


M 


Fig.  30.  Isoplanatic  Angle  Measurements,  Lick  Observatory 
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Fig.   31.   Spatial  Coherence  Leng-th  Measurements, 

Lick  Observatory 
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Fig.   32.   Echosounder  Trace,  Lick  Observatory 
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Fig.   33.   Echosounder  Trace,  Lick  Observatory 
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Fig.   34.   Echosounder  Trace,  Lick  Observatory 
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The  strong  correlation  between  the  echosounder  data  around 
0400  to  0500  local  standard  time  and  the  spatial  coherence 
length  ( ro )  measurements  at  1200  universal  time,  combined 
with  the  lack  of  correlation  between  the  isoplanatic  angle 
( 8o )  and  the  spatial  coherence  length  ( ro )  measurements 
indicate  that  the  lower  atmospheric  and  surface  turbulent 
layers  are  dominating  the  atmospheric  profile  during  this 
time  period. 

Data  collected  on  10  April  1987,  again  illustrate  the 
strong  correlation  between  the  three  atmospheric 
measurements  made.  A  comparison  of  the  isoplanatic  angle 
(So )  and  the  spatial  coherence  length  ( ro )  measurements 
during  the  time  interval  of  0700  and  1300  universal  time 
indicate  steady  turbulent  conditions  in  the  upper 
atmosphere  and  greatly  varying  turbulent  conditions  at 
lower  atmospheric  levels.  This  is  evident  in  Figures  35 
and  36  by  the  consistent  values  of  the  isoplanatic  angle 
( 9o )  during  the  time  period  compared  with  the  steady 
increase  and  eventual  decline  of  the  spatial  coherence 
length  ( ro )  values  during  the  same  time  interval.  The 
trace  variations  in  Figure  36  during  the  hours  of  0800  and 
1200  universal  time  are  indicative  of  a  period  of 
diminishing  lower  atmospheric  or  surface  turbulence 
followed  by  the  onset  of  an  increasingly  turbulent  period 
around  1200.  This  turbulent  trend  is  distinctively 
supported  by  the  echosounder  data  in  Figures  37  through  44. 
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Fig.   35.   Isoplanatic  Angle  Measurements,  Lick  Observatory 
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Fig.   36.   Spatial  Coherence  Length  Measurements, 

Lick  Observatory 
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Fig.   37.   Echosounder  Trace,  Lick  Observatory 
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Fig.   38.   Echosounder  Trace,  Lick  Observatory 
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Fig.   39.   Echosounder  Trace,  Lick  Observatory 
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Fig.   40.   Echosounder  Trace,  Lick  Observatory 
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Fig.   41.   Echosounder  Trace,  Lick  Observatory 
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Fig.   42.   Echosounder  Trace,  Lick  Observatory 
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Fig.   43.   Echosounder  Trace,  Lick  Observatory 
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Fig.   44.   Echosounder  Trace,  Lick  Observatory 
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The  large  200  mm  spatial  coherence  lengths  around  1100 
universal  time  are  consistent  with  the  low  turbulence 
evident  in  the  echosounder  profiles  around  0300  local 
standard  time. 

Only  acoustic  echosounder  measurements  were  made  on  the 
upper  roof  of  Spanagel  Hall  at  the  Naval  Postgraduate 
School.  Data  runs  on  26  and  27  April  1987,  include  both 
echosounder  data  and  the  associated  Ct 2  plots.  Of 
particular  interest  is  the  plot  of  11:30  on  the  27th  of 
April  in  Figure  45,  which  shows  the  maritime  boundary 
inversion  layer  around  100  m  being  perturbed  by  convective 
plumes  at  lower  altitudes.  Sixteen  additional  echosounder 
output  traces  were  included  in  Appendix  D  to  exhibit 
typical  atmospheric  activity.  Many  of  the  plots,  such  as 
the  1 4: 45,  15:30,  and  the  17:30  of  26  April  and  the  10:00, 
11:15.  13:00,  14:15,  15:15,  16:15,  and  the  17:15  of  27 
April,  have  convective  plumes  which  are  prevalent  whenever 
a  heat  flux  between  the  surface  and  atmosphere  exists.  The 
plots  of  18:15  and  19:15  on  the  26th  of  April  and  of  13:00, 
18:45,  and  20:00  on  the  27th  of  April  clearly  show  the 
passage  of  the  neutral  event  which  is  encountered  when  the 
atmospheric  and  surface  temperature  difference  becomes 
negligible.  Finally,  the  plot  of  16:30  on  26  April  can  be 
associated  with  strong  winds  which  exhibit  a  somewhat 
uniform  return  for  all  altitudes  across  the  entire  15 
minute  interval  of  the  trace. 
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Fig.   45.   Echosounder  Trace  and  Ct 2  Plot,  NPS  Roof 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A  high  frequency,  5kHz,  phased  array  echosounder 
constructed  to  measure  low  level  turbulence  appears  to  work 
within  the  10  to  135  meter  altitude  range.  When  using  the 
device,  detailed  profiles  of  the  short  range  atmospheric 
density  fluctuations  were  obtained.  The  profiles  were 
found  to  correlate  very  well  with  the  measurements  of  the 
isoplanatic  angle  (9o )  and  the  spatial  coherence  length 
( ro )  during  periods  of  simultaneous  operation.  This  short 
range  echosounder,  when  used  in  conjunction  with  the  other 
atmospheric  measuring  devices,  is  an  invaluable  tool.  The 
graphical  output  provides  us  with  a  more  complete 
description  of  the  atmosphere  that  can  be  used  to  calculate 
On2 ,  the  atmospheric  index  of  refraction  structure 
parameter . 

Further  software  developments  such  as  the  incorporation 
of  Doppler  and  Fast  Fourier  Transform  routines  into  the 
echosounder  program  will  determine  the  radial  velocity 
profile  of  the  return  signal.  Such  an  improvement  will 
determine  the  velocity  of  the  probed  air  masses  and  plot 
these  echosounder  traces  as  a  function  of  color  intensity. 

Other  areas  of  further  reasearch  include  the  design  and 
testing  of  different  array  patterns.  One  such  pattern,  the 
hexagonal   array,   is   already   undergoing   tests.    As  the 
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speaker  array  evolves,  it  is  also  evident  that  the 
enclosure  must  follow  suit  to  accomodate  the  new  array 
design  and  associated  beam  patterns.  Further  improvements 
in  software  routines  are  also  inevitable.  It  is  always 
desirable  to  store  on  a  disc  all  data  collected  at  a  site. 
Presently,  the  floppy  disc  capacity  is  inadequate  for  the 
storage  of  more  than  two  hours  of  data.  Furthermore,  if 
FFT  and  Doppler  routines  are  to  be  added,  the  computational 
speed  of  the  computer  will  inhibit  the  pulse  repetition 
rate.  A  special  data  acquisition  technique  called  "Direct 
Memory  Access"  may  then  have  to  be  introduced  into  the  HP 
computer  to  allow  simultaneous  data  collection  and 
processing . 

Finally,  this  product  may  be  used  in  various 
applications  not  already  addressed.  One  such  use  may  be  to 
measure  windshear  at  airports.  Other  researchers  have 
expressed  interest  in  using  the  echosounder  to  measure 
arctic  atmospheric  conditions  during  meteorological 
surveys.  Assuredly,  as  the  device  evolves  into  a  more 
compact  and  highly  efficient  instrument,  its  range  of 
application  will  continue  to  expand. 
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(ANTENNA  BEAM  PATTERN  FOR  AN  N  ELEMENT  LINEAR  ARRAY 

GCLEAR 

Ri=.0381 

Fr=5000 

Ka=(Ri*2*PI*Fr)/343 

INPUT  "ENTER  THE  «  OF  DRIVERS  "  ,N 

D=.0762 

Kd=(2*PI*D*Fr)/343 

Kd2=Kd/2. 

Num_decades=4 

Scale=1 0ANum_decades 

Thetastep=.001 

GCLEAR 

VIEWPORT  20,80,20,80 

WINDOW  -20,105,-2,1  10  . 

AXES  100/(2*Nun_decades ) , 1 00/( 2 . *Num_decades ) ,0,0,20,20 

MOVE  10,103 

LABEL  "PATTERN  FOR  ";N;"  ELEMENTS" 

CSIZE  3 

FOR  1=1  TO  5 

MOVE  -15,(25»(I-1  ))-1 
M=(  10-((6-I  )*10)) 
LABEL  M 
NEXT  I 
LDIR  PI/2 
MOVE  -15,40 
LABEL  "DECIBELS" 
LDIR  0 
FOR  Theta=Thetastep  TO  PI/2  STEP  Thetastep 

Arg=Ka*SIN(Theta) 

Arg2=2*FNJKArg) 

Arg2=Arg2/Arg 

Kds=Kd2*SIN(Theta) 

Htheta=SIN(Kds*N )/< N*SIN(  Kds ) ) 

Htheta=Htheta*Arg2 

Htheta=Scale*Htheta*Htheta 

Htheta=LGT( Htheta )/Num_decades 

X=Htheta*SIN(Theta ) 

Y=Htheta*COS(Theta> 

PLOT  X* 1 00  ,Y* 1 00 

(PRINT  Theta.Htheta 
NEXT  Theta 
END 
I 
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450    DEF  FNJKX  ) 

460    !  J1  BESSEL  FUNCTION  10  MAY  ,1984;  DLW 

470    IF  X>3.  THEN  GOTO  Xgt3 

480    X3=X/3. 

490    X3=X3*X3 

500    Fx=.S+X3*(-.56249985+X3*( .  2 1 093573+X3*( - . 03954289+X3* 

(  .  004433 1 9+X3*<-. 00031 76 1+X3*. 00001 109) ) ) ) ) 
510    J1=Fx*X 
520    RETURN  J1 
530  Xgt3:   X3=3./X 
540   Fx=.79788456+X3*( . 000001 56+X3*( .  01 659667+X3*( .0001 71 05+X3* 

(-.0024951 1+X3#(  .001 13653-X3*.00020033))  ))) 
550    Tx=X-2.35619449+X3*( . 1 24996 1 2+X3*( .0000565+X3*( -.00637879+ 

X3*( .00074348+X3*( . 00079824-X3* . 000291 66 ) ) ) ) ) 
560    J1=Fx*COS(Tx )/SQR(X ) 
570    RETURN  J1 
580    FNEND 
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APPENDIX  B 
LCDR  BUTLER'S  BEAM  PATTERN  PROGRAM 
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PRINT  "**  BEAM  PATTERN  COMPUTATION  AND  PLOT  PROGRAM  **' 
PRINT  "***********************************************•' 
t 
t 

PRINT  "*»***  SETJJP  *****" 

i 

i 

PRINT  CHR$  (10) 

PRINT  "THIS  PROGRAM  CONTROLS  THE  FOLLOWING  APPARATUS" 
PRINT  "TO  MEASURE,  DISPLAY,  PLOT  AND  STORE  THE" 
PRINT  "DIRECTIONAL  CHARACTERISTICS  OF  A  SINGLET," 
PRINT  "DOUBLET,  TRIPLET,  AND  QUADRUPLET  ARRAY." 
t 

PRINT  "  INSTRUMENT  ADDRESS 

PRINT  "  

PRINT  "  HP_3478  MULTIMETER  (VOLTS) 

PRINT  "  HP  3478  MULTIMETER  (RES) 


724 
723 
709 


PRINT  "  DATA  AQUISITION/CONTROL 
t 

DISP  "INITIALIZE  INSTRUMENTS" 
f 

REM  ***  PRESET  AND  INITIALIZE  INSTRUMENTS*** 

OUTPUT  709  ;"RS"  !  RESET 

OUTPUT  723  ;"F3R5T1"  !  FREE  RUN  30K0HMS 

OUTPUT  724  ;"F2RAT1"  !  FREE  RUN  AUTORANGE  AC  VOLTS 

f 

DISP  "INITIALIZE  PROGRAM" 
t 

OPTION  BASE  1 

COM   ANGLEH50)  ,CALC_ANGLE(  150)  ,NAME$[801 

DEG 

PRINTER  IS  701 

INTEGER  I.CL.N 

DIM  R(360) ,VOLT(360) ,OHMS_DATA(360 ) ,X(360) ,Y(360) fK(360) 

MVOLT=0 

MASS  STORAGE  IS  ".DATA1" 

DISP  "NUMBER  OF  DATA  POINTS  DESIRED?" 

BEEP  10,300 

INPUT  N 

DISP  "ENTER  DISC  DATA  FILE  NAME" 

INPUT  NAMES 

DISP  "WANT  NEW  DATA  <Y/N>?" 
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430  BEEP  10,300 

440  INPUT  0$ 

450  IF  Q$="Y"  THEN  GOTO  480 

460  IF  Q$="N"  THEN  GOTO  920 

470  GOTO  430 

480  REM  ***  DEFAULT  PARAMETERS  *** 

490  ! 

500  !  MUST  INSERT  FOR  LINES  390,  400,  410  AND  810 

510  !  PRIOR  TO  START  OF  RUN 

520  ! 

530  CREATE  NAMES, 360, 16 

540  ASSIGN*  1  TO  NAMES 

550  NAME$="RUN_1" 

560  FULLSCALE=500  !  LOCK_IN  MVAC 

570  CONV=360/10000  !  CONVERT  OHMS  TO  DEGREES 

580  SETJJP:  !  ***  SETJJP  MEASUREMENT  PARAMETERS  *** 

590  PRINT  CHR$  (10) 

600  PRINT  "1. GRAPH  TITLE :" ;NAME$ 

610  REM  ***  MODIFY  IN  PRECEEDING  LINE  NHEN  NECESSARY  *** 

620  DISP  "MANUALLY  ENTER  DESIRED  FREQ  FOR  SPEAKER  DRIVER" 

630  DISP  "ENTER  0  IF  NO  CHANGE" 

640  DISP  "ENTER  1  TO  CHANGE  TITLE" 

650  BEEP  10.600 

660  INPUT  C  !  CHOICE 

670  IF  C  =  0  THEN  GOTO  690 

630  IF  C=1  THEN  GOTO  830 

690  PRINT  CHR$  <10> 

700  OUTPUT  709  ;"CLS0"  !  ROTATE  CW 

710  FOR  1=1  TO  N 

720  OUTPUT  709  ;"CLS1"  !  POWER  TO  MOTOR  ON 

730  WAIT  .5  !  MILLISECONDS 

740  OUTPUT  709  ;"OPN1"  !  POWER  TO  MOTOR  OFF 

750  WAIT  4000 

760  ENTER  723  ;  OHMS_DATA(I) 

770  ENTER  724  ;  VOLT(I) 

780  PRINT  "I  IS  NOW"; I 

790  PRINT  "OHMS_DATA=";OHMS_DATA(D 

800  PRINT  "VOLT=";VOLT<I> 

810  NEXT  I 

820  FUR  1=1  TO  N 

830  PRINT*  1  ;  OHMS_DATA< I) , VOLT( I) 

840  NEXT  I 

850  ASSIGN*  1  TO  * 

360  GOSUB  PLOT_ROUTINE . 

370  GOTO  910 

380  ENTER  NAMES  !  NEW  TITLES 

890  GOTO  600 

900  PRINT  "DONE" 
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910  END 

920  DISP  "WHAT  IS  DATA  RUN  DESIRED'?" 
930  BEEP  10,500 
940  INPUT  NAMES 
950  ASSIGN*  1  TO  NAMES 
960  FOR  1=1  TO  N 

970  READ*  1  ;  OHMS_DATA< I) , VOLT< I) 
930  NEXT  I 
990  ASSIGN*  1  TO  * 
1000  GOTO  360 
1010  ! 
1020  ! 

1030  PLOT_ROUTINE:  !  ***  PLOTS  POLAR  DATA  *** 
1040  ! 
1050  ! 

1060  PLOTTER  IS  705 
1070  GRAPH 
1080  GCLEAR 
1090  PEN  1 

1100  LIMIT  50,210,20,180 
1110  SCALE  -60,60,-60,60 
1120  XAXIS  0,10,-40,40 
1130  YAXIS  0,10,-40,40 
1140  CSIZE  3 
1150  LINE  TYPE  1 
,1160  FOR  L=0  TO  40  STEP  10 
1170  MOVE  0,L 
1180  LORG  6 
1190  LABEL  L-40 
1200  NEXT  L 
1210  LINE  TYPE  1 
1220  CSIZE  5 
1230  MOVE  0,50 
1240  ! 

1250  !  INSERT  LABELS  1250  AND  1280  PRIOR  TO  START  OF  RUN 
1260  ! 

1270  LABEL  ;"5  X  5  ACOUSTIC  ARRAY  BEAM  PATTERN" 
1230  CSIZE  3 
1290  MOVE  0,-45 
1300  LABEL  ;"R=3.60  M'V'MIKE  VT=" , VOLT< I ) , "INPUT=2 .5  V",' 

","FREQ=5.0  KHZ" 
1310  CSIZE  4 
1320  MOVE  45,0 
1330  LORG  5 
1340  LABEL  ;"MIKE" 
1350  LINE  TYPE  4 
1360  FOR  L=10  TO  40  STEP  10 
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1370  MOVE  L.O 

1380  FOR  1=0  TO  360  STEP  5 

1390  DRAW  L*COS  (I)  ,L*SIN  (I) 

MOO  NEXT  I 

1410  NEXT  L 

1420  LINE  TYPE  1 

1430  PEN  UP 

1440  PEN  2 

1450  FOR  1=1  TO  N  !  BEGINS  BEAM  PATTERN  COMPUTATION 

1460  IF  VOLT<I)>MVOLT  THEN  MVOLT=VOLT< I > 

1470  NEXT  I 

1430  FOR  1=1  TO  N 

1490  K=0HMS_DATA(I)/10000*360 

1500  K1=K 

1510  IF  OHMS_DATA(I)<10  THEN  K=K1+360/295 

1520  IF  OHMS_DATA(I)>10500  THEN  GOTO  1610 

1530  VECTOR=20*LGT  (VOLT( I)/MVOLT)+40 

1540  LINE  TYPE  1 
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0  IF  VECT0R<1  THEN  VECTOR=0 


1560  PRINT  I, K, VECTOR 

1570  X(I)=VECTOR*COS  (K) 

1580  Y(I)=VECTOR*SIN  <K> 

1590  PLOT  X(I),Y<I) 

1610  NEXT  I 

1615  PEN  UP 

1620  RETURN 


71 


APPENDIX  C 
ECHOSOUNDER  SOFTWARE 

The  driving  force  of  this  system  is  found  within  the 
HP200  Series  Computer.  The  software  and  computer  system 
were  responsible  for  controlling  and  monitoring  every 
operational  phase  of  the  hardware  components.  The  program 
which  accomplishes  this  task  was  entitled  "ACRDR"  and  is 
listed  in  Reference  3.  This  program  was  modeled  after  a 
program  written  by  Walters  [Ref.  20],  but  each  program 
performs  a  distinctly  different  computational  task.  This 
program  was  written  in  HP  Basic  4.0  and  was  compiled  by  an 
Infotek  BC203  Basic  Compiler  to  enhance  the  speed  of 
execution . 

The  program  "ACRDR"  is  easily  broken  down  into  a  number 
of  blocks  and  subroutines  which  perform  specific 
operations.  These  sections  are  outlined  in  a  flowchart 
(Fig.  46).  Each  block  is  straightforward  in  its  purpose, 
and  the  program  is  designed  to  be  as  helpful  to  the  user  as 
possible.  As  a  prologue  to  the  actual  program  code,  there 
is  a  listing  of  all  the  program  variables  with  a  short 
desription  of  their  use.  Such  a  listing  familiarizes  a 
user  with  the  computations  to  be  made  and  also  provides 
quick  references  for  any  future  modifications  to  be 
Implemented . 
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The  initialization  procedures  which  set  all  the 
parameters  necessary  for  data  collection  follow  the 
variable  definition  section.  These  parameters  control  such 
things  as  the  contrast  between  background  and  return  signal 
and  the  setting  of  the  computer's  internal  clock.  All 
options  of  the  HP3314A  Function  Generator  are  decided  by 
the  user  at  this  phase.  Another  initialization  involves 
preparing  the  Infotek  AD200  Analog  to  Digital  Converter  for 
operation.  Finally,  internal  arrays  are  dimensioned, 
integer  variables  are  defined  and  the  computer  function 
keys  are  redefined  to  suit  the  "ACRDR"  program. 


Initialization 
Procedures  for: 

HP  331 4A 

AD  200  A/D  Converter 

Internal  Clock 

Constants 


I 


Computation  of 
Parameters 

Speed  of  Sound 

Offset 

Range 

Pulse  Length 

Attenuation  Coefficient 


I 


Prepare  Screen  for 
Data  Display 


Trigger  Pulse  from 
HP  331 4A 


Send  Cy   Plot  to  the 
Printer 


Compute  Values  of  CT 
and  Plot 


Send  Plot  to  the 
Printer 


Collect  Return  and 
Perform  Reduction 
Routines,  Store  Data  if 
Required 


Compute  the  Average 
Noise  Level 


Fig.   46.   Flowchart  of  Fundamental  Program  Operations 
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Now  that  the  system  is  ready  to  operate,  certain 
constants  must  be  determined  based  upon  atmospheric 
conditions  at  the  time  of  data  collection.  From  the 
ambient  temperature,  a  speed  of  sound  in  air  is  determined 
by  the  equation  in  Kinsler,  et  al ,  [Ref.  14:  p.  106]. 

C  =  Co  [  1  +  (T/273)  ]i/2 

where 

C  is  the  local  speed  of  sound, 

Co  is  the  speed  of  sound  in  air  at  0°  Celcius,  (331.6 
m/s ) ,  and 

T  is  the  local  temperature  in  degrees  Celsius. 
Using  this  value  and  a  user  decided  sampling   frequency  and 
data   point   total,   the   maximum   range   of   the  traces  is 
determined  by  the  relation, 

RANGE  =  C/2  •  l/( SAMPLE  FREQUENCY)  •  N  , 

where 

RANGE  is  the  range  of  the  echosounder  in  meters, 
SAMPLE  FREQUENCY  is  the  user  desired  rate  of  data 

acquisition,  and 

N  is  the  number  of  data  points  per  trace. 

Typically,  the  values  of  many   of   the   parameters   are  not 

varied.    The  value  of  N  above  is  almost  always  taken  to  be 
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the  allowed  maximum  of  16301  data  points  and  the  speed  of 
sound  can  be  safely  estimated  to  be  340  m/s  at  room 
temperature.  By  defining  a  typical  sampling  frequency  of 
20000  Hz,  a  range  of  approximately  135  meters  is  obtained. 

Another  parameter  determined  is  the  atmospheric 
attenuation  coefficient.  This  calculation  is  done  in  a 
subroutine  obtained  from  Reference  21  and  is  necessary  in 
this  program  for  the  calculation  of  the  temperature 
structure  parameter,  Ct 2 .  Other  computations  made  prior  to 
program  execution  include  the  pulse  length,  wavenumber,  and 
the  D-C  offset  of  the  equipment. 

At  this  point,  the  computer  is  finally  ready  for  data 
collection.  The  screen  setup  displays  a  distance  versus 
time  plot  and  the  internal  clock  of  the  computer  is 
synchronized  with  the  pulse  from  the  HP3314A  Function 
Generator.  Following  transmission  of  the  pulse,  the  return 
signal  is  received,  digitized,  stored  in  an  array  and  data 
reduction  commences.  A  block  averaging  technique  is  used 
in  which  a  block  of  data  points  the  size  of  the  number  of 
emitted  cycles  is  summed  and  averaged.  Returns  are  plotted 
as  darkened  areas  on  the  echosounder  traces  with  the 
intensity  of  the  darkened  area  directly  related  to  the 
magnitude  of  the  return  signal.  This  same  procedure  is 
used  in  the  computation  of  Ct 2  with  some  very  important 
differences.  Aside  from  simply  correcting  for  n:r2 
spherical  divergence,  Ct 2  is  corrected  for  electronic  gain, 
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the  ratio  of  power  returned  to  power  transmitted,  the 
efficiency  of  the  speakers  as  transmitters  and  receivers, 
the  area  of  the  speaker  array,  the  atmospheric  attenuation, 
the  pulse  length,  the  temperature,  the  scattering  cross 
section  per  unit  volume  at  a  specific  range  and  frequency 
and  finally  the  effective  aperture  factor  of  the  antenna. 
Additionally,  Ct 2  (Range)  is  averaged  for  a  particular 
altitude  over  15  minute  intervals.  Finally,  after  each 
pulse  is  reduced,  a  corresponding  mean  square  noise  level 
is  determined.  This  is  done  by  averaging  each  block 
average  at  maximum  range  until  at  least  ten  values  have 
been  used  in  the  average.  An  upper  limit  of  five  over  the 
average  noise  level  (this  corresponds  to  voltage 
fluctuations  on  the  order  of  ID-7  volts)  is  set  on  the 
rountine  to  avoid  averaging  any  strong  return  signals  or 
anomalies  such  as  passing  aircraft.  After  the  ten  values 
are  averaged,  every  subsequent  pulse  is  averaged  into  all 
the  preceding  noise  levels  and  removed  from  each  subsequent 
return  signal. 

Ultimately,  after  each  15  minute  interval,  or  at  the 
user's  request,  the  display  terminal  image  is  printed.  At 
this  time  the  Ct 2  computations  are  conducted  and  plotted  on 
the  screen  and  printer  as  range  versus  logi o  of  Ct 2 .  Upon 
completion  of  the  printing  of  these  plots,  the  system 
begins  the  data  processing  for  the  next  15  minute  interval. 
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There  are  certain  options  built  into  the  program  to 
allow  a  user  to  change  various  aspects  of  operation.  The 
function  keys  allow  the  user  to  change  the  local 
temperature,  sample  frequency  or  intensity  factor  during 
program  execution.  Additionally,  the  user  can  quit  or 
restart  the  program,  print  the  partial  trace  on  the  screen 
or  elect  to  save  a  trace  on  a  floppy  disc.  The  save 
routine  is  invoked  for  the  subsequent  15  minute  interval 
after  the  appropriate  function  key  is  depressed.  Saving  a 
future  trace  may  be  ackward,  especially  if  a  user  would 
like  to  keep  an  interesting  trace  which  is  presently  on  the 
terminal.  This  problem  cannot  be  readily  solved  unless 
each  trace  is  recorded  to  disc  without  user  intervention. 
At  present,  this  is  not  done  because  only  eight  traces  (two 
hours  of  data)  can  be  written  to  a  floppy  disc  before  it  is 
full. 
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APPENDIX  D 
ECHOSOUNDER  OUTPUT 
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